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ABSTRACT 

We present Submillimeter Array observations of the z=3.91 gravitationally lensed broad absorption 
line quasar APM 08279+5255 which spatially resolve the 1.0 mm (200 /im rest-frame) dust continuum 
emission. At 0'.'4 resolution, the emission is separated into two components, a stronger, extended one 
to the northeast (46 ± 5 mJy) and a weaker, compact one to the southwest (15 ± 2 mJy). We have 
carried out simulations of the gravitational lensing effect responsible for the two submm components 
in order to constrain the intrinsic size of the submm continuum emission. Using an elliptical lens 
potential, the best fit lensing model yields an intrinsic (projected) diameter of ~80 pc, which is not 
as compact as the optical/near-infrared (NIR) emission and agrees with previous size estimates of the 
gas and dust emission in APM 08279+5255. Based on our estimate, we favor a scenario in which the 
200 /im (rest-frame) emission originates from a warm dust component (Td=150-220 K) that is mainly 
heated by the AGN rather than by a starburst (SB). The flux is boosted by a factor of ~90 in our 
model, consistent with recent estimates for APM 08279+5255. 

Subject headings: galaxies: high-redshift — galaxiesigravitationally lensed — submillimeter igalaxies — 
galaxies: individual: (APM 08279+5255) 



1. INTRODUCTION 

APM 08279+5255 (=APM08) is a_strongly lensed 
broad absorption line (BAL) quasar (Irwin et al .11 19981 : 
[Lewis et al. 199 8) with a very po werful active galactic 
nucleus (AGN; ISoifer et al.l Il987l ). This combination 
makes it a n extremely br ight object despite its redshift 
of 2^=3.91 (|Downes et al.| [l999). Its bolometric luminos- 
ity is thought to be ~5xl0^ L©, amplified by up to 
a factor of 100 by a foreground galaxy which has yet 
to be identified. Due to the large amplification, broad 
absorption lines hav e been detected even in the X-ray 
([Chartas et al.|[2Q02l ). Ground-based, Chandra and HST 
observations have suggested that APM08 c onsists of at 
least two compon ents separate d by ^ 0^^4. Ilbata et all 
(jl999l . =199) and lEgami et all (|2000l . =E00) later de- 
tected a third image (their image C), and argued that it 
is likely a third lensed ima g:e of the bac kground QSO in- 
stead of the lensing galaxv. lLewis et al] (|200 2b) obtained 
optical spectra of the three images using HST/STIS, and 
showed that the three spectra are quite similar. This in- 
dicates that the image C is indeed a third lensed image 
of APM08. Due to the magnification by the lens, sev- 
eral molecular lines have been detecte d in APM 08 (e.g., 
Weifi et a l. 2007; Guehn et al. 2007; Riechers et al. 2006; 
Wagg et a l. 2006; Garcia-Buri llo et all 2006: Wagg et al. 

2OO5I : iLewis et al 2002a: Downes et al 1999). suggest- 
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The 

dust m ass is estima ted to be Mdust=5xl0^ M© 
([Downes et al.l 119991 : IWfi et al.l 120071 . W07). However, 
despite the increasing amount of mm-data, no tight con- 
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straints have been set on the size of the dust and/or 
gas emission which may help to characterize the heating 
source of the dust and gas, i.e., whether it is in the form 
of an AGN and/or starburst (SB). Only recently have 
W07 presented some indirect evidence through bright- 
ness temperature arguments that the molecular line and 
dust emission originate from a compact region with a ra- 
dius of 100-200 pc. In this Letter, we present observations 
from the Submillimeter Array (SMA) with sufficient reso- 
lution (0'.'4) to resolve the 1.0 mm dust continuum emis- 
sion in APM08, and lensing models that constrain its 
intrinsic source size. 

2. OBSERVATIONS 

We observed the 1.0 mm (200 fim rest-frame) con- 
tinuum emission in APM08 with the SMA using seven 
of the eight 6 m antennas in two configurations, with 
baselines ranging from 20 to 500 meters, in November 
and December 2006. This results in a synthesized beam 
of 0'.'42x0'.'35 at a position angle (PA) of 141° using 
robust weighting. The LO frequency was 302.4 GHz 
(1.0 mm), with each sideband centered 5 GHz to either 
side. This frequency setup was chosen so that the [N U] 
line (1461 GHz rest frame) falls within the LSB. However, 
no significant singal of this line has been found above a 
3(7 threshold of 9 Jy beam~^ km s~^, integrated over 
a velocity range of 1000 km s~^ around the expected 
[N n] frequency. We used Titan as a flux calibrator, 
3C111, 3C454.3, J0927+390 and 3C279 as bandpass cali- 
brators, and J0927+390 as a gain calibrator. J0753+538 
has been frequently observed throughout the night as 
a test source to assess the quality of the gain calibra- 
tion. The phase rms on the gain calibrator is <40 de- 
grees, resulting in a positional uncertainty of <0'.'l in the 
test map of J0753+ 538, which is sim ilar to theoretical 
noise considerations (|Reid et al.lll988l ), and to a 'seeing' 
of <0.03'' (=difference between the apparent size of the 
test quasar and the beam size). The accuracy of the 
flux calibration is around 20%. Weather conditions were 
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good with an atmospheric opacity of r(225 GHz) =0.06- 
0.12 (recorded at the nearby Caltech Submihimeter Ob- 
servatory). The data were written into the MIR format 
originahy developed at the Owens Vahey Radio Obser- 
vatory (Scovihe e t al.l [199 3). and reduced using SMA- 
specific MIR tasks in IDL. Fohowing cahbration, the data 
were exported to FITS and imaging and further analysis 
were done using GILD AS. Both 2 GHz sidebands were 
merged together in the uv-plane to image the 1.0 mm 
continuum emission, resulting in an effective rms noise 
of 1.7 mJy beam~^. 

3. RESULTS 

The 1.0 mm continuum in APM08 is clearly detected 
at high angular resolution (~0'.'4) and high sensitivity 
(SNR>9; Fig.[T]). Using robust weighting, the two lensed 
images, labeled NE (north-eastern) and SW (south- 
western) respectively (Fig[T]), are clearly separated (by 
~0'.'4) from each other. The image positions listed in Ta- 
ble [T] are measured after the image was deconvolved with 
the synthesized beam. Both absolute positions and the 
separation are consistent with previous results within the 
uncertainties (199, EOO). A comparison of the robustly 
and naturally weighted (angular resolution ~0'.'6) maps 
indicates similar total fluxes of ~60±12 mJy (Table [1]). 
Taking the (single dish) fluxes of 24 mJy at 1.35 mm and 
of 75 mJy at 0.850 mm a nd assuming a spectral index of 
2.46 (iLewis et al.l (|1998f )). a flux of 50 mJy is expected 
at 1 mm which is, within the uncertainties, in agreement 
with the SMA value. This indicates that no flux has 
been resolved out in our data and the dust emission at 
200 /im rest-frame cannot be very extended, similar to 
recent results in APM08 by W07. 

The deconvolved size (=FWHM) of the emission in NE 
appears to be slightly more extended (by ^30%) than 
the synthesized beam (Table [T]). An extension in NE is 
further supported in comparing the peak flux with the 
spatially integrated flux (Table [T]). While they differ by 
a factor of ~2 for NE, they are very similar in SW. This 
suggests that SW is still unresolved at the high angu- 
lar resolution of the robustly weighted map. We derive 
a peak flux ratio of (NE/SW)peak ^1.4 and a spatially 
integrated flux ratio of (NE/SW)int -3. The latter is 
significantly higher than the value found in the (com- 
pact) optical/NIR emission of ~1.2 (e.g., 199, EOO), indi- 
cating differential lensing effects. The morphology of the 
200 /im rest-frame dust emission is in good agreement 
with that of the CO(l-O) emission of APM08 recently 
published bv lWalter et al. I (1200 6). Based on their C0(1- 
0) map fFig. 4 in lWalter et al. h 2006). the NE/SW peak 
flux ratio is estimated to be ^1.5 and the FWHM of NE 
is estimated to be 0^3-0^4, which is slightly larger than 
their beamsize, while SW is clearly unresolved. 

4. SIMULATIONS OF THE LENSING EFFECT 

Previous lensing models for APM08, obtained by fit- 
ting the (compact) NIR emission in EOO, fail to reproduce 
the (extended) emission seen in our submm observations, 
as in their extended emission models the southern image 
turns out to be too strong and too extended. There- 
fore, we have run new simulations^ of the gravitational 

^ Assumed cosmology: r2M=0.3, I^r — 0.0 and Qy^O.? fo r a Hub- 
ble constant of Ho=71km s~^ Mpc~^ (jSoergel et al.ll2003l V giving 
an angular-size distance of 1)^4=1.47 Gpc and 1^'=7 kpc. 



l ensing effect usin g the code developed by Krips & Neri 
(|Krips et al.ll2005L available in GILDAS). We concentrate 
our simulations on an elliptical lens potential (=EP) sim- 
ilar to previous studies (e.g., EOO), although we also tried 
a single isothermal sphere (SIS) as a lens potential. Al- 
though the SIS model fits our submm data, we discarded 
the SIS models based on the assumption that APM08 is 
a three-image system and because of the inability of the 
SIS model to produce an odd number of lensed images. 

We varied the parameters of the lensing potential as 
well as the size of the unlensed source in a reasonable pa- 
rameter space (see Table [2j) until the best agreement be- 
tween the observed and simulated maps was found, i.e., a 
model with the lowest reduced test (<1). We used the 
position of the two lensed images (see Table [T]), their peak 
flux ratio (NE/SW=1.4±0.2) as weh as their (lensed) 
shape (FWHMc^0'.'5±0'.'l in NE and FWHM<0'.'2 in SW) 
as constraints. We augmented this set of constraints with 
those from the (compact) optical/NIR emission (e.g., 199, 
EOO; Fig. [Tj), accounting also for a third lensed image 
close to NE. This results in 9 observational constraints 
versus 8 free lensing parameters. We processed the simu- 
lated lensed images through a spatial filter defined by the 
uv-coverage of the observations. This is a more robust 
approach than simply convolving the simulated data with 
the synthesized beam of the observations as it also ac- 
counts for possible resolution effects (see also iKrips et al.l 
2005). However, based on this set of observational con- 
straints, we could not find any lensing solution that re- 
produces the optical/NIR and our submm observations. 
We doubt that this inconsistency is due to mismatched 
angular resolutions between the optical/NIR and our ob- 
servations. When starting with an elliptical potential 
that yields three lensed images and fits the optical/NIR 
data, NE always turns out to be either too strong or not 
extended enough compared to SW to match the submm 
data, or, equivalently, SW is either too extended or not 
strong enough compared to NE to be consistent with the 
submm emission. Also, models that fit the submm emis- 
sion always make the lensed optical/NIR image NE too 
strong or too extended compared to SW. 

A solution to reproduce the optical/NIR and submm 
data with the same lensing model is to allow different 
positions between the unlensed optical/NIR and submm 
sources of emission. This actually increases the number 
of free lensing parameters to 10. However, we have to 
handle the observed positions of the lensed images in the 
optical/NIR and submm case as independent parameters 
as well, so that the observational constraints also increase 
to 11. Even though the positional accuracy of the submm 
and the optical/NIR observations do not yet allow us to 
substantiate the positional offset (~0'.'08±0'.'07) between 
the (lensed) optical/NIR and submm emission and thus 
also between the unlensed positions, this is still a reason- 
able approach since it yields a lens model that fits both 
optical/NIR and submm data. This lens model, called 
EP(3), is plotted in Fig. [2^ andEJo and listed in Table El 
For comparison. Fig. [2t and [2]i respectively show sim- 
ulations in which the unlensed submm source is smaller 
and larger than the best-fit size. The lens parameters 
and the position of the unlensed submm source are the 
same as in Fig [2^ and [2)3. As the optical/NIR source 
is very compact in its intrinsic size, the simulation in 
Fig. [2]3 is the one to be compared with the optical/NIR 
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observations. The NE/SW ratio for this compact case, 
however, disagrees with the optical/NIR observations, 
and the peak separation is clearly wider than what is 
observed. To correctly reproduce the optical/NIR data, 
the position of the unlensed source has to be shifted by 
-0.003'' (=21 pc; Fig.Et). Such an offset between opti- 
cal/NIR and submm emission could be a consequence of 
asymmetrically distributed submm emission in APM08. 
The (best-fit) intrinsic source diameter of the submm 
emission is —0.012" (=84 pc) and the magnification fac- 
tor —90, consistent with recent results from WOT. 

5. NATURE OF THE DUST EMISSION 

WOT and B06 (jBeelen et al.ll2006[ ) both present a two 
temperature model for the dust emission in APM08 with 
an extended 'cold' {Td=^0-70 K) and a compact 'warm' 
(T(i= 160-220 K) component. The warm dust component 
is assumed to be heated by the AGN up to a maximum 
radius of 60-130 pc (WOT). W07 also estimate that the 
dust heating from the AGN can still reach —65 K at a 
radius of 350 pc and may dominate the cold component 
as well. However, this estimate has to be considered as 
an upper limit because of potential self-screening effects 
in the AGN. Therefore, a significant contribution to the 
cold component from more extended star formation can- 
not be completely excluded. 

By looking at the dust SED in W07^, the contribution 
of the 'warm' component to the 200/im rest-frame emis- 
sion is —70%. The diameter estimate of 80 pc based 
on our data lies within the maximum possible diam- 
eter of up to 120-260 pc in W07 for the 'warm' dust 
component. Even when accounting for the apparent off- 
set between AGN and dust emission, this may indicate 
that the 200 /im rest-frame emission in APM08 is domi- 
nated by the warm dust component although a fractional 
contribution of the cold component cannot be entirely 
excluded. Following previous assumptions (e.g., W07) 
and a similar argumen tation as used for Arp 220 (see 
iDownes fc Eckartll2007l ). a SB within a diameter of 80 pc 
seems unlikely also in APM08 leaving the AGN as the 
main suspect for the dust heating at 200 /im: similar 
to Arp 220, APM08 has a very high intrinsic luminos- 
ity (for m=90) of -2x1O^2L0 at 200 fim (W07) and, 
hence, an emission surface brightness of >1O^^L0 kpc~^, 
which is 60 times the luminosity of the nearby SB galaxy 
M82 in a -4000 times smaller volume (M82's SB disk 
has a radius of —650 pc; iGarcia-Burillo et al.l [2000). 
Although single super-star clusters (SSCs) can reach 
such a high surface brightness of ^10^^L(t^ k pc~^ (e.g. 
iMartm- Hernandes. Sch aerer. fc S auvagd l2005l ). APM08 
would need to harbor >1000 SSCs that would fill al- 
most its entire central region (typical SSC diameters are 
—1 pc). This seems to be a rather unlikely scenario. 
Even if we consider the uncertainties of our models, we 
can clearly constrain the intrinsic radius to be ^150 pc 
(see Fig.[2li), unless the dust 'disk' in APM08 had a very 

^ As WOT uses a grey body fit to the dust SED and B06 only a 
black body, we base our discussion on WOT. 



small inclination angle (z<30°), which is, however, not 
supported by CO observations of WOT who favor i^70°. 
Also in this more conservative case, the inferred surface 
brightness of APM08 still exceeds that of any SBs seen in 
the present-day universe by more than an order of mag- 
nitude, likely discarding the SB scenario for APM08. 

Based on the flux and the deconvolved size determined 
from the submm-observations, the instrinsic brightness 
temperature at rest of APM08 can be estimated to be 
Tb^30 K. A comparison to the dust temperature of 
220 K of the warm component suggests a dust opac- 
ity of r^O.lb (with r=— ln(l-Tb/Tdust))- Assuming the 
dust temperature of the cold component of 65 K yields a 
more conservative upper limit of the opacity of <0.6 at 
200 fim rest-frame, translating to an even smaller upper 
limit when converted to 1.3 mm. This upper limit is sig- 
nificantly lower than that found for the AGN i n Arp 220 
of >0.7 at 1.3 mm bv lDownes fc EckartI (|2QQ7r ): the dust 
emission in Arp 220 has been suggested to be very similar 
to that of APM08. This may indicate that the dust emis- 
sion in APM08 could be either less dense and/or more 
clumpy than around the AGN of Arp 220. However, 
given the similar (intrinsic) dust masses of ^lO^M© and 
the similar intrinsic diameters of ~80 pc in both sources, 
a lower density in APM08 may seem less likely. 

6. SUMMARY &; CONCLUSION 

We have detected the 200 jam rest-frame continuum 
emission in APM08 using the SMA with an angular 
resolution of ~0'.'4. The two (main) lensed images are 
clearly separated with a combined flux of ~60±12 mJy. 
Simulations of the gravitational lensing effect in this sys- 
tem yield a diameter of the intrinsic (submm) continuum 
emitting region of ~80 pc and magnification factor of 
90. Our data and simulations seem to be only consistent 
with a lensing scenario including a third lensed image 
of APM08 close to NE if the position of the compact 
optical/NIR emission and the position of the extended 
submm emission are offset from each other before lens- 
ing by ~0'.'003 (=21 pc). Further (sub)mm observations 
may be beneficial to determine whether such a positional 
offset is indeed necessary or more complicated lens po- 
tentials have to be considered. Given our size estimate, 
we favor a scenario in which the 200 /im emission origi- 
nates from a warm dust component that is supposed to 
be mainly heated by the AGN rather than by a SB. 
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TABLE 1 

Observational Parameters of APM 08279+5255 



Component 


n 


n 


qpeak J, 
^Imm 

(mJy beam"-*^) 


cint c 
^Imm 

(mJy) 


Deconvolved Size 
maj. xmin.,P.A. 

(//^Z/^O) 


Total^ 


-0.10+0.03 


-0.10+0.03 


34+2 


60+12 


(0.8x0.7)+0.1,(52+5) 




0.00+0.04 


0.00+0.04 


21+2 


46+5 


(0.5x0.2)+0.1, (50+20) 




-0.28+0.06 


-0.26+0.06 


15+2 


15+2 


<0.2 



^ The offsets are with respect to Q;j2ooo=08h31m41.69s and 5j2ooo=52d45ml7.5s. The positional errors comprise the statistical errors 
from the Gaussian fit to the data in the uv-plane and uncertainties from the calibration. ^ Peak flux. Errors are based on the statistical 
noise in the map.^ Spatially integrated flux. The flux errors include the statistical uncertainties from the flt and in the case of the total 
flux also the calibrational uncertainties, which are estimated to be ~20%.^ Based on the naturally weighted map and an elliptical gaussian 
flt. ^ Based on the robustly weighted map and an elliptical (circular) gaussian flt for NE (SW). 



TABLE 2 

Parameters of the best-fit gravitational lens model. 









PA^ 












n 




n 


n 




n 




EP(3)" 


0.31 


0.008 


-30 


0.21 


(-0.218,-0.158) 


0.012 


90 



^ Elliptical potential with three lensed images; we allowed the unlensed source position to be variable between the optical/NIR 
and our submm data (see text). ^ Einstein radius ^e; tested parameter range: ^e=0'.'1-0'.'4^ Ellipticity e; tested param- 
eter range: e=0-0.2^ Positional angle PA (from N to E): tested range —90 to +90°^ Core radius 6c'-> tested parameter 
range: ^c=0-0.3''^ The position of the unlensed source is with respect to the lensed component NE. The lens position is at 
(AQ;,A(5)=(-0''210,-0''l50).g (9s=intrinsic source diameter; tested range: O'.'OOl-O'.'S (=0.007-2.1 kpc)^ m=magnification factor 



6 



Krips et al. 




0.5 -0.5 -1 0.5 -0.5 -1 

RA - offset (") RA - offset (") 



Fig. 1. Left: HST image of APM08 observed with the WFPC2/F814W filter (STScI archive). The black contours represent the 1mm 
continuum emission. Right: Continuum image at the observed wavelength of 1.0 mm of APM08 with robust weighting. Contours start at 
±3(7 in steps of 2cr; la corresponds to 1.7 mJy beam~^. The yellow crosses mark the (deconvolved) position of the two lensed images; they 
are thus slightly more separated than the peaks of the contour image. The image is referenced to component NE (Table 
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Fig. 2. — Simulations (color scale and white contours) based on an elliptical lens model with three lensed images: EP(3). The simulated 
emission is overlaid with the observed emission (black contours). Contours are the same as in Fig. ^ The model with the best-fit source 
size for the submm data are shown in 2a) (unfiltered) and in 2b) (spatially filtered). 2c) and 2d) columns show the best-fit model with the 
smallest and largest unlensed source sizes, while 2e) shows the model for the optical/NIR case (smoothed to 0(^15). The models in 2c) and 
2d) refers to the (unlensed) submm position; the (unlensed) optical/NIR position is offset by ~ 0(^003 in 2e) from the (unlensed) submm 
position (see Section |4]). 



